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Abstract 

We present an exact solution for electronic Raman scattering in a single-band, strongly correlated material, includ- 
ing nonresonant, resonant and mixed contributions. Results are derived for the spinless Falicov-Kimball model, 
employing dynamical mean field theory; this system can be tuned through a Mott metal-insulator transition. 
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Electronic Raman scattering is an important probe 
of electronic excitations in materials. It has been used 
to examine different kinds of charge and spin excita- 
tions in a variety of different materials, ranging from 
Kondo insulators [1,2], to high temperature supercon- 
ductors [3,4], to collossal magnetoresistance materi- 
als [5] . Inelastic light scattering involves contributions 
from scattering processes that depend on the incident 
photon frequency (so-called mixed and resonant con- 
tributions) and processes that are independent of the 
incident photon frequency (so-called nonresonant con- 
tributions). There has been much theoretical work on 
this problem. In the strong-coupling regime, a pertur- 
bative approach has been used, and has illustrated a 
number of important features of resonant scattering 
processes [6,7]. The nonresonant case has also been ex- 
amined, and an exact solution for correlated systems 
(in large spatial dimensions) is available for both the 
Falicov-Kimball [8] and Hubbard [9] models. Here we 
concentrate on an exact solution of the full problem for 
the Falicov-Kimball model including all resonant and 
mixed effects. 



The Falicov-Kimball model Hamiltonian [10] is (at 
half filling) 
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and includes two kinds of particles: conduction elec- 
trons, which are mobile, and localized electrons which 
are immobile. Here c\ (a) creates (destroys) a con- 
duction electron at site i, Wi is the localized electron 
number at site i, U is the on-site Coulomb interaction 
between the electrons, and t* is the hopping integral 
(which we use as our energy unit). The symbol d is the 
spatial dimension, and {ij) denotes a sum over all near- 
est neighbor pairs (we work on a hypercubic lattice). 
The model is exactly solvable with dynamical mean 
field theory when d — » oo [11,12] (see [13] for a review). 

Shastry and Shraiman [6] derived an explicit formula 
for inelastic light scattering that involves the matrix 
elements of the electronic vector potential for light with 
the many-body states of the correlated system. The 
expression for the Raman response is 
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for the scattering of electrons by optical photons (the 
repeated indices a and /3 are summed over). Here e^f) 
refer to the initial (final) eigenstates describing the 
"electronic matter" , Z is the partition function, and 
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is the scattering operator constructed by the current 
j a = J2u 9s(k) /dk a c^Ck, and stress-tensor 7^,3 = 
J2k <9 2 e(k) / dk a dkpc^Ck operators, with e(k) the band 
structure and Ck the destruction operator for an elec- 
tron with momentum k. 

Instead of calculating the matrix elements in Eq. (3) , 
we use a diagrammatic technique to calculate all con- 
tributions to the Raman response function x(^)i which 
is defined by 
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Our calculations include effects from nonresonant di- 
agrams, from resonant diagrams, and from so-called 
mixed diagrams. We have evaluated these expressions 
for the Stokes Raman response, with an incident pho- 
ton frequency u>i, an outgoing photon frequency u) , 
and a transfered photon frequency 57 = Wj — 0J o - The 
procedure is complicated, and involves first computing 
the response functions on the imaginary time axis, then 
Fourier transforming to imaginary frequencies, and fi- 
nally performing an analytic continuation to the real 
axis [14]. 

We analyze three different symmetries for the inci- 
dent and outgoing light. The Ai g symmetry has the 
full symmetry of the lattice and is measured by tak- 
ing the initial and final polarizations to be e l = e° = 
(1,1,1,...) (we assume nearest-neighbor hopping only). 
The Big symmetry is a d-wave-like symmetry that in- 
volves crossed polarizers: e l = (1,1,1,...) and e° = 
(—1,1,-1,1,...). Finally, the i?2g symmetry is another 
d-wave symmetry rotated by 45 degrees; with e 1 = 
(1,0,1,0,...) and e° = (0,1,0,1,...). 

The total Raman response function is the sum of the 
nonresonant, mixed, and resonant contributions and 
has a complicated form. It turns out that the ^4i g sector 
has contributions from nonresonant, mixed, and reso- 
nant Raman scattering, the Bi g sector has contribu- 
tions from nonresonant and resonant Raman scatter- 
ing only, and the B2 g sector is purely resonant [8]. It 
is educational to consider the contributions of the bare 
diagrams, which can be summed up and rewritten in 
the following form: 
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Ak(u>) — -ilmGk(w — iS), Gk(co) is the momen- 
tum-dependent single-electron Green's function, and 
f(u>) — 1/[1 + exp(Pco)] is the Fermi distribution func- 
tion. 

In general, the bare response function in Eq. (5) is 
a function of the frequency shift = uji — ui , of the 
incoming photon frequency u>j and the outgoing pho- 
ton frequency ui ; it can be enhanced when one or both 
of the denominators of the Green's functions of the 
effective vertex (expression on the last line) are reso- 
nant. In the latter case, we have a so-called "double" 
or "multiple resonance" [15]. The full response func- 
tion also includes vertex renormalizations. But the to- 
tal (reducible) charge vertex for the Falicov-Kimball 
model does not diverge, and hence it does not introduce 
any additional energy denominators or "resonances" ; 
it only leads to a renormalization of the total Raman 
response. 

The Falicov-Kimball model on a d — 00 hypercu- 
bic lattice has a Mott transition into a pseudogap-like 
phase at half filling when U = a/2- We examine the 
system on the insulating side of the Mott transition at 
[7 = 3. 

Our results for the total Raman response as a func- 
tion of the transfered frequency Q for different tem- 
peratures appear in Fig. 1 for the incident photon fre- 
quencies u>i = 2.5 and 00, and in Fig. 2 for u)i = 2. 

For the u>j = 00 case (thin lines in Fig. 1) we have a 
pure nonresonant response that is nonzero only in the 
Ai g and Bi g channels. One can see, that in both chan- 
nels all lines, which correspond to different tempera- 
tures, cross at a characteristic frequency fl ~ U/2 (isos- 
bestic point) where the Raman response is independent 
of temperature. For the Falicov-Kimball model at half- 
filling, all of the temperature dependence in Eq. (5) en- 
ters from the Fermi distribution function f(u>). In the 
insulating phase, the rapidly varying parts of the Fermi 
distribution functions are located in the gap regions of 
the two-particle density of states, and the temperature 
dependence is strongly reduced when these gaps be- 
come symmetric for some values of the transfered and 
incident frequencies. 

In Ref. [8], an isosbestic behavior was observed for 
the nonresonant response only in the Bi g channel, 
but in Fig. 1 it is also seen in the Ai g channel when 
the response is plotted on a logarithmic scale. When 
the incident photon frequency decreases (thick lines 
in Fig. 1), we also have a non-zero Raman response in 
the i?2g channel and the shape of the Raman response 
is changed. In particular, a sharp peak appears at the 
double resonance located at Q = w», and the full re- 



2 



Q[t] 



S 1x10"" K' 

a 1 



CO 



1x1 0'- 

1x10 s 

1x10"' 

1x1 s 

1x1 O c 
1x10" 
1x10" 
1x1 0" ! 
1x10"" 
. 1x1 0" E 
1x10 s 
1x10"' 
1x1 0" f 
1x1 0" E 



s 



I 
I 

I 

I 



T=1 
T=0.5 
T=0.2 
T=0.05 



0,5 



!2[t] 





I 


1 I 1 I 


1 






























y 




















- 
























i; 
n 










il 


T=1 








/; 


T=0.5 










T=0.2 




/ 






T=0.05 




r 












-t 


i . i 


■ 





Q[t] 

Fig. 1. Isosbestic behavior of the resonant Raman response for 
U = 3 and uJi — 2.5 (thick lines) and uJi — oo (nonrcsonant 
response, thin lines). Different lines correspond to different 
temperatures T = 1, 0.5, 0.2, 0.05. Note that the bJi = 2.5 
curves all cross at two isosbestic points: one close to U /2 and 
another close to LJi. 



sponse is not just an enhancement of the nonresonant 
features (which are apparent when the incident pho- 
ton frequency becomes large), but the shape of the 
response can change dramatically due to resonant ef- 
fects. This is most apparent when the incident photon 
energy is close to U. Also, the initial UJi — oo isosbestic 
point is shifted and a second isosbestic point appears 
in the double resonance area of Q w Wj. With a fur- 
ther decrease of the incident photon frequency, both 




Fig. 2. The same plot as in Fig. 1, but for the case uJi — 2. 

isosbestic points approach one another (Fig. 2) and 
then disappear when Ui ~ U/2. 

In Fig. 3, we plot the total Raman profile at various 
transfered frequencies 51 as a function of the incident 
photon frequency u>i. Note that when Q, is larger than 
the energy of the charge-transfer excitation (Q. > 17) 
we only observe the double resonance at Ui = Q and 
there are no additional features in the resonant profile. 
When Q decreases and moves into the charge-transfer 
peak region, a resonant enhancement of the charge- 
transfer peak at u>i w U is observed. Its location and 
width change with a further decrease of the transfered 
frequency and it almost disappears when Q lies between 
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Fig. 3. Raman response at various values of the transfcrcd 
frequency H in steps of 0.2 as a function of the incident photon 
frequency uJi for T — 0.5. All curves start at oji — Q. 

the charge-transfer and low-energy peaks, and is then 
restored when f2 moves into the low-energy peak re- 
gion, where it has a double-peak structure for the Bi g 
and £>2 g channels. So, we observe a joint resonance of 
the charge-transfer and low-energy peaks, and this res- 
onance of a low-energy feature due to a higher-energy 
photon has been seen in the Raman scattering of some 
strongly-correlated materials. 

In conclusion, we have performed an exact calcula- 
tion of the electronic Raman response function for a 
strongly correlated system in the insulating phase and 
predict four interesting resonant features: (1) the ap- 



pearance of a double-resonance peak, (2) the nonuni- 
form enhancement of non-resonant features due to res- 
onance, (3) the appearance of two isosbestic points in 
all channels, and (4) a joint resonance of the charge 
transfer and low energy peaks when the incident pho- 
ton frequency is on the order of U. It will be interest- 
ing to see whether these features can be seen in future 
experiments on correlated systems. 

Acknowledgments: We acknowledge support from 
the U.S. Civilian Research and Development Foun- 
dation (CRDF Grant No. UP2-2436-LV-02). J.K.F. 
also acknowledges support from the National Science 
Foundation under Grant No. DMR-0210717. T.P.D. 
acknowledges the NSERC, PREA and the Alexander 
von Humboldt foundation for support. 



References 

[I] P. Nyhus, S. L. Cooper, and Z. Fisk, Phys. Rev. B 51, 
15626 (1995). 

[2] P. Nyhus, et al., Phys. Rev. B 52, R14308 (1995); Phys. 
Rev. B 55, 12488 (1997). 

[3] K. B. Lyons, et al. Phys. Rev. Lett. 60, 732 (1988): S. 
Sugai, S. Shamoto, and M. Sato, Phys. Rev. B 38, 6436 

(1988) ; P. E. Sulewsky, et al. Phys. Rev. B 41, 225 (1990); 
R. Liu, et al. J. Phys. Chem. Solids 54, 1347 (1993). 

[4] G. Blumberg, et al. Phys. Rev. B 53, 11930 (1996). 

[5] H. L. Liu, et al. Phys. Rev. B 60, R6980 (1999). 

[6] B. S. Shastry and B. I. Shraiman, Phys. Rev. Lett. 65, 
1068 (1990); Int. J. Mod. Phys. B 5, 365 (1991). 

[7] A. V. Chubukov and D. M. Frenkel, Phys. Rev. B 52, 9760 
(1995); Phys. Rev. Lett. 74, 3057 (1995); D. K. Morr and 
A. V. Chubukov, Phys Rev. B 56, 9134 (1997). 

[8] J. K. Freericks and T. P. Devereaux, Condens. Matter 
Phys. 4, 149 (2001); Phys. Rev. B 64, 125110 (2001). 

[9] J. K. Freericks, T. P. Devereaux, and R. Bulla, Phys. Rev. 
B 64, 233114 (2001); J. K. Freericks, et al. Phys. Rev. B 
67, 155102 (2003). 

[10] L. M. Falicov and J. C. Kimball, Phys. Rev. Lett. 22 
(1969) 997. 

[II] W. Metzncr and D. Vollhardt, Phys. Rev. Lett. 62 (1989) 
324. 

[12] U. Brandt and C. Mielsch, Z. Phys. B— Condens. Mat. 75 

(1989) 365; 79 (1990) 295; 82, (1991) 37. 

[13] J. K. Freericks and V. Zlatifi, Rev. Mod. Phys. 75 (2003) 
1333. 

[14] A. M. Shvaika, O. Vorobyov, J. K. Freericks, and 
T. P. Devereaux, unpublished. 

[15] R. M. Martin and L. M. Falicov, in Light Scattering in 
Solids, edited by M. Cardona (Springer- Verlag, New York, 
1975). 



4 



